Draft version September 15, 2010 

Preprint typeset using LSTgX style emulateapj v. 1 1/10/09 



CONTINUED COOLING OF THE CRUST IN THE NEUTRON STAR LOW-MASS X-RAY BINARY KS 1731-260 
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ABSTRACT 

Some neutron star low-mass X-ray binaries (LMXBs) have very long outbursts (lasting several years) which 
can generate a significant amount of heat in the neutron star crust. After the system has returned to quiescence, 
the crust then thermally relaxes. This provides a rare opportunity to study the thermal properties of neutron star 
crusts, putting constraints on the thermal conductivity and hence the structure and composition of the crust. 
KS 1731-260 is one of only four systems where this crustal cooling has been observed. Here, we present a 
new Chandra observation of this source approximately 8 years after the end of the last outburst, and 4 years 
since the last observation. We find that the source has continued to cool, with the cooling curve displaying a 
simple power-law decay. This suggests that the crust has not fully thermally relaxed yet, and may continue to 
cool further. A simple power law decay is in contrast to theoretical cooling models of the crust, which predict 
that the crust should now have cooled to the same temperature as the neutron star core. 
Subject headings: stars: neutron — X-rays: binaries — X-rays: individual: KS 1731-260 



1. INTRODUCTION 

Neutron stars in transient LMXBs provide a rare observa- 
tional opportunity to study the thermal properties of the neu- 
tron star crust. In some of these transient systems accretion 
outbursts last many years, as opposed to the more typical 
weeks to months. These quasi-persistent systems are particu- 
larly interesting as during the long outbursts the crust should 
be heated out of ther mal equilibrium with the rest of the star 
dRufledge et alJ2002l) . Thus, once the source returns to quies- 
cence, the crust thermally relaxes. The resulting cooling curve 
depends on a number of key properties of th e crust includ- 
ing the composition and structure of the crust dRutledge et alj 
l2002tlShternin et al]2007HBrown & Cumming|2009l) . and the 
crust thi ckness (which is dependent on the mass and ra dius of 
the star. lLattimer et alJll994UBrown & Cumming||2009l) . 

Such crustal cooling has now be en observed in four 
neutr o n star transients: KS 1731-260 (I Wijnands et al.l 120011 
2001 ICackett et all l2006h MXB 1659-29 dWiinands et alj 
20031 12004 ICackett et all 120061 l2008h . EXO 0748-676 
dDegenaar et all 120091 120101) and XTE J1701-462 
jFridriksson et alj(2010I) . KS 1731-260 was in outburst for 



12.5 years, returning to quiescence in 2001 (Wijnands et al 
120011). For a detailed history of this source, see Cackett et al 
(2006). It was quickly realized that such a long outburst 
should allow for observable crustal cooling, hence, a moni- 
toring campaign using Chandra and XMM-Newt on followed 
dWiinands et alJlIoMIIOOl ICackett et alj|2006l) . A cooling 
cu rve covering the first ~ 4 years of quiescence was presented 
by ICackett et all d2006l) . The source flux initially dropped 
rapidly: a factor of ^2.5 during the first yea r in quiescence 
and a factor of ~ 5 in th e first 1000 days (Wijnan ds et alj 
2002| ICackett et all 120061) . The cooling curve was well 
fit by either an exponential decay to a constant level or a 
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simple power-law decay dCackett et alj 120061 [2008). The 
rapid rate of cooling and the inferred cold temperature of the 
neutron star core, led to conclusions that the crust must have 
a high thermal c onductivity and requ i red enhanced levels 
of core cooling dWiinands et alj 120021: ICackett et all 120061) 
based on comparing with the theoretical cooling models of 
iRufledge et alj d2002l) . Motivated by these cooling curves, 
theoretical c rust cooling models for KS 1731-260were 
calculated by IShternin et al.l (120071) . These authors rule out 
a low thermal conductivit y for the crust, but d o not r equire 
enhanced core cooling. iBrown & Cu mming (2009!) also 
calculated thermal relaxation models for neutron star crusts, 
finding a high thermal conductivity for the crust, and a low 
impurity parameter (i.e., the dispersion in the charge of the 
ions in the crust is low). 

Here, we present a new Chandra observation of 
KS 1731-260, performed ~ 8 years after the end of the out- 
burst, and 4 years since the previous Chandra observation. 
The data are consistent with further cooling of the neutron 
star, continuing along a power-law decay. 

2. ANALYSIS AND RESULTS 

The newest Chandra observation of KS 173 1-260 was per- 
formed in two separate pointings. One 31 ks segment was 
performed on 2009 May 17 (ObsID: 10037) and the other 28 
ks segment was performed on 2009 May 19 (ObsID: 1091 1). 
The source was at the default aimpoint on the ACIS-S3 chip, 
which was operated in FAINT mode. Given several recent 
changes to the calibration of Chandra, we have opted to also 
reanalyze all the previous Chandra observations with the lat- 
est calibration and software. Moreover, we reanalyze the 
XMM-Newton observations with the latest calibration files and 
softwar e. For full details of the previous observations, please 
refer to ICackett et ail d2006l) . 

2. 1 . Chandra data reduction 
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The Chandra data were all analyzed using CIAO (v 4.2) and 
CALDB (v 4.2.2). We used a circular source extraction re- 
gion with a radius of 3", and the background extraction region 
was an annulus with inner radius of 7"and outer radius 25". 
The psextract tool was used to extract the spectra, and 
mkacisrmf and mkarf were used to create the response 
matrices. 

2.2. XMM-Newton data reduction 

The XMM-Newton data were analyzed using the XMM 
Science Analysis Software (v 9.0.0). The observation data 
files were reprocessed using the emproc and epproc tasks. 
To check for high levels of background flaring we extracted 
lightcurves for all events with > 10 keV and pattern = for 
the MOS, and 10 - 12 keV and pattern = for the PN. In 
all observations there was some significant background flar- 
ing. We excluded all times where the > 10 keV lightcurve 
was greater than 2 c/s for the MOS and 4 c/s for the PN. We 
extracted spectra using evselect with a circular source ex- 
traction region of radius 10", and a circular background ex- 
traction region of radius 1' taken from a source free region 
close to the source. For the MOS we filtered with patterns - 
12, and for the PN, we used patterns 0-4 and flag = 0. Re- 
sponse matrices were generated using rmf gen and arf gen. 

2.3. Spectral analysis 

We fit the spectra using XSPE C (ver. 12, lArnaudlll996l) . 
following similar p rocedures to ICackett et al.l (120061) and 
ICackett et al.l (12008) . The spectra were modeled with an ab- 
sorbed neutron star atmosphere model. We used the phabs 
model for Galactic ab sorption, and the n sa model for the neu- 
tron star atmosphere (Za vlin et al.| [l996). Throughout, we fix 
the neutron star radius at 10 km, and the neutron star mass 
at 1.4 M©. The normalization of the nsa model is given by 
l/D 2 , where D is the dista nce to the source in pc. Here, we 
assume a distance of 7 kpc (M uno et al.12 000). i.e. normaliza- 
tion =2.041 x 10~ 8 p c~ 2 . As we showed for MXB 1659-29 in 
ICackett et al.l (12008). the distance assumed only shifts the fit- 
ted temperatures up or down, and does not affect the cooling 
timescales derived. Here, if we assume a distance of 5 kpc, 
the temperatures are all about 10% lower than for D = 7 kpc, 
where as if we assume D = 9 kpc then the temperatures are 
approximately 7% higher than for D = 7 kpc. 

We fit all the spectra simultaneously with the absorption 
column density, A^h, the same for all spectra, and a free pa- 
rameter in the fit. The neutron star atmosphere effective tem- 
perature is allowed to vary between epochs. For the first 
XMM-Newton observations (013795201/013795301) the two 
sets of spectra had their parameters tied between the spectra. 
Similarly, for the last Chandra observation (10037/1091 1) we 
also tie the parameters between the two spectra. Given the 
low count rates there are not enough counts per bin to use \ 2 
statistics when fitting, and therefore we use the W-statistic in 
XSPEC to fit the unbinned spectra. 

The results of the spectral fitting are given in Table Q] and 
the evolution of the effective temperature is shown in FigureQ] 
All uncertainties quote d and plotted are at th e la level of con- 
fidence. Note that in Cacke tt et alj (1200 6) the uncertainties 
quoted were at the 90% level of confidence, not 1 a as stated 
in the text. The results clearly show a decrease in the neu- 
tron star atmosphere temperature over time, with the newest 
observation the coldest yet. 
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FIG. 1 . — Effective temperature for an observer at infinity for KS 1731—260 
over approximately 3000 days from the end of the last outburst. The solid 
line shows a power-law fit to these temperatures, while the dashed line 
shows the best-fitting exponential decay to a constant level. The dotted 
line shows the best- fitting model using the crustal cooling simulations from 
Brow n & Cumminel <2009l) . 

2.4. Count rate analysis 

One potential problem with our interpretation of crustal 
cooling is that there could be a power-law spectral compo- 
nent that cannot be detected here because of the low num- 
ber of counts in the spectra. Non-thermal power-law compo- 
nents are common in many quiescent neutron star spectra (e.g. 
Uonker et al.|[2004b . Any power-law component could poten- 
tially change the shape of the cooling curve of the thermal 
component. It should be noted that the highest quality spec- 
tra we have of KS 173 1-260 is the first Chan dra observation. 
This observation was discussed in detail by IWijnands et al. 
( 120011) . These authors show that the spectrum is consistent 
with a thermal spectrum only, and a power-law is not required 
statistically. However, if they add a power-law to their model 
they find that it only contributes ^15% to the flux. Due to the 
low count rate of the later observations, here, we use count ra- 
tios between different bands to assess whether there could be 
a significant power-law component present. One can only do 
this with observations from the same telescope (Chandra and 
XMM-Newton have different effective areas), thus we analyze 
just the raw counts from the six Chandra observations. 

Figure [2] shows Chandra count rate versus time. This 
lightcurve also shows a power-law decay like the tempera- 
tures from spectral fitting. We find that only the very first 
observation has a significant detection in the 3-10 keV band 
where the power-law component would dominate. To study 
the spectral evolution of the source, we therefore define a 
hardness color ratio between the count rate in the 1.5 - 3.0 
keV band and the 0.5 - 1.5 keV band. In Figure[3]we show a 
hardness-intensity diagram that displays the evolution of this 
color with the full 0.5 - 10 keV count rate. We also show 
the predicted color and count rate evolution for three different 
spectral models: (i) a neutron star atmosphere, (ii) a neutron 
star atmosphere plus power-law model where the power-law 
always contributes 20% to the unabsorbed 0.5 - 10 keV flux 
(we assume a spectral index, T = 2), (iii) a neutron star at- 
mosphere plus power-law (V = 2) model where the power-law 
flux always contributes 40% of the unabsorbed 0.5 - 10 keV 
flux (note that all three models assume Galactic absorption of 
Nh = 1.3 x 10 22 cm" 2 ). The data points are most consistent 
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TABLE 1 

Neutron star atmosphere spectral fitting parameters 





2428 
(CXO) 


013795201/301 
(XMM) 


3796 
(CXO) 


3797 
(CXO) 


0202680101 
(XMM) 


6279 
(CXO) 


5468 
(CXO) 


10037/10911 
(CXO) 


MJD 

N H (10 22 cnT 2 ) 


51995.1 


52165.7 


52681.6 


52859.5 
1.30 ±0.06 


53430.5 


53500.4 


53525.4 


54969.7 



kT°° (eV) 103.2 ± 1.7 88.9± 1.3 75.5±2.2 73.3±2.3 71.0± 1.8 66.0±4.5 70.3 ±2.1 63.1 ±2.1 

fboi a (10~ 13 ergcirT 2 s" 1 ) 4.2±0.3 2.3±0.1 1.2±0.1 1.1 ±0.1 0.95±0.09 0.72±0.21 0.92±0.11 0.60±0.08 



NOTE. — All uncertainties are lcr. The nsa atmosphere model is used and the mass and radius are fixed in all fits to 1.4 Mq and 10 km respectively. kT^ is the effective temperature 
for an observer at infinity, i.e. the gravitationally redshifted effective temperature. The normalization of the atmosphere was set for a distance = 7 kpc. ObsIDs of each observation are 
indicated at the top, with CXO (Chandra) and XMM {XMM-Newton) denoting the observatory. 
a Bolometric flux calculated from the model over the 0.01 - 100 keV range 



> 0.0010 



100 1000 
Time since end of outburst (doys) 

FIG. 2. — 0.5-10 keV count rates for all Chandra observations of 
KS 1731—260. The solid line shows a power-law fit to these count rates. 

with a simple cooling neutron star atmosphere, but, model (ii) 
where there is a small contribution from a power-law cannot 
be ruled out. Any power-law contribution increases the color 
(for a given count rate) compared to the color from a neutron 
star atmosphere only. A harder power-law spectrum (with a 
slope closer to 1 than 2) would also lead to an increased color. 
Thus, any contributions from a power-law component must 
be at a low level (less than a few tens of percent). Therefore, 
there must still be significant neutron star cooling. Also note 
that this count rate analysis does not take into account any 
change in the effective area of the detector over time. There 
has been a known increase in contaminant on the Chandra 
ACIS detector since the mission launch, leadin g to a decrease 
in the sensitivity at lowest energies over time (Marshall et al. 
2004). This would artificially harden the color used here, thus 
the colors shown here should be taken as maximum values. 

To further test the affect of a constant fraction of power-law 
flux, we fit the spectra with an absorbed neutron star atmo- 
phere plus power-law model. We fix the power-law with slope 
r = 2 and the normalization for each observation is set to give 
a 0.5 - 10 keV unabsorbed flux that is 20% of the value found 
from fitting a neutron star atmosphere alone. The resulting 
temperatures still follow a simple power-law decay cooling 
curve, with the same slope as before. This can be easily un- 
derstood, as reducing the thermal flux by 20% will simply 
just reduce the measured temperature by approximately 5% 
as F oc T 4 . 

2.5. Cooling curves 

The latest Chandra observation analyzed here was taken 
approximately 4 years after the nearest observations, and ap- 
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FIG. 3. — Hardness-intensity diagram for Chandra observations of 
KS 1731-260. The hardness color is the ratio of the 1.5 -3.0 keV count 
rate to the 0.5 - 1.5 keV count rate. This is plotted against the count rate for 
the 0.5 - 10 keV band. The solid line shows the predicted count rates for a 
cooling neutron star atmosphere. The dashed line shows the predicted count 
rates for a neutron star atmoshere and power-law spectrum, where the power- 
law (r = 2) always contributes 20% to the unabsorbed 0.5 - 10 keV flux. The 
dotted line shows a neutron star atmosphere plus power-law spectrum (r = 2) 
where the power-law flux always contributes 40% to the unabsorbed 0.5 - 10 
keV flux. 



pears to show continued cooling of the neutron star. We com- 
pare the three observations taken around MJD 53500 with this 
newest observation. These three observations (the fifth, sixth 
and seventh quiescent observations) were all taken within 100 
days of each other, and all have consistent temperatures at the 
lcr level. A weighted average of the temperatures from these 
three observations gives kT^ = 70.3 ±1.3 eV. Compared to 
the latest Chandra observation, where kT^ = 63.1 ±2.1 eV, 
the two temperatures differ at the 3a level of confidence, sug- 
gesting continued cooling of KS 1731-260. 

Before this latest observatio n, the cooling curve could be 
well fit by either a power-law (Cackett et al. 200 8|) or a n ex- 
ponential decay to a constant level dCackett e t al. 2006J). We 
tested both forms here. We fitted a power-law of the form 
y(f) = a(f-fo)' 3 to the effective temperatures^, where we chose 
fo to be midday on the las t day the source wa s observed to be 
active (MJD 51930.5, see lCackett et al.l2006l) . The best fitting 
parameters are a = 174.7 ±1.3 eV and /3 = -0.125 ± 0.007, 
giving xt = 0.33, P x = 0.92. This best-fitting power-law is 
shown in Figure [T] as a solid line, and fits the cooling curve 
well. For comparison, we also test an exponential decay of 

2 Note that as F oc 7"* the flux decay curve can simply be calculated from 
the fit to the temperatures 
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the form y(t) = aexp[-(t-to)/b] + c. The best fitting parame- 
ters are a = 39.8 ± 2.3 eV, b = 418 ±70 days and c = 67.7 ± 1 .3 
eV, giving xt = 2.0, P x = 0.04. The simple power-law decay 
is therefore a better fit to the cooling curve, further indicating 
continued cooling of KS 1731-260. 

We also fitted numerica l crust cooling simulations from 
IBrown & Cummins! (2009) to the data. These calculations 
by lBrown& Cumming (2009) suggest that the lightcurve of 
a cooling crust is expected to be a broken power-law which 
flattens to a constant at late times, set by the temperature of 
the neutron star core. The initial power-law decay is set by 
the temperature profile in the outer crust. The power-law 
break occurs when there is a transition in the solid from a 
classical to quantum crystal at a depth close to neutron drip. 
The time of this power-law break is set by the thermal diffu- 
sion time to the depth of this transition. In these models the 
three variable parameters are (i) the crust impurity parameter 

which is given by Qi mp = »^^~]»,(Z, — (Z)) 2 and measures 

i 

the dispersion in the charge of the nuclides (Z) in the crust, 
(ii) the core temperature, T c and (iii) the temperature at the 
top of the crust, T b . The best-fitting model (dotted line in 
Fig [TJ has Q imp = 4.0, T c = 5.4 x 10 7 K, and T h = 2 x 10 8 K 
(at a column of 1 x 10 12 g cm" 2 during the outburst). This 
model giv es y 2 = 3.3. These values are quite similar to the 
best-fit in Bro wn & Cumm ing (2009) who found Qi mp = 1.5, 

T c = 4.6 x 10 7 K, and T b = 2.5 x 10 8 K. Thus, the core temper- 
ature here is slightly colder and the impurity parameter in the 
crust is slightly higher than the previous best-fitting values. 



3. DISCUSSION AND CONCLUSIONS 

We have presented a new Chandra observation of the neu- 
tron star low-mas X-ray binary KS 1731-260 in quiescence. 
This observation extends the monitoring of this source to ap- 
proximately 8 years after the end of the most recent outburst. 
The observation suggests that the neutron star has continued 
to cool with an effective temperature (for an observer at in- 
finity) of kT^ = 63.1 ±2.1 eV compared to the preceding 
epoch where kT^ = 70.3 ±1.3 eV. Previously, the cooling 
curve could be well fit by both a po wer-law decay or an expo - 
nential decay to a constant level (ICackett et al.ll2006L 12008). 
Here, we find that the latest observation remains on the sim- 
ple power-law decay with a slope of -0.125 ±0.007, indicat- 
ing that the crust may be continuing to cool and has not yet 
reached thermal e quilibrium with the core. As has been dis- 
cussed previously (IWiinands et al.ll2002l: I Cackett et al1l2006t 
IShternin et al.l l2007t IBrown & Cumming! 120091) . the cooling 
of KS 1731-260 indicates a crust with high thermal conduc- 
tivity. This implies that the crust has a low impurity parame- 
ter, in other wor ds, the dispersion in the ch arge of the ions in 
the crust is low (Br own & Cum ming 2009, find Q imp < 10 is 
a robust upper limit). 

One potential problem is that the spectrum may not be 
purely thermal, but that there could also be some contribution 
from non-thermal power-law emission that cannot be detected 
due to low number of counts in the spectra. Recent observa- 
tions of another cooling neutron star XTE J1701-462 have 
shown an increase in flux for approximately 100 days on top 



of a cooli ng curve, likely due to i ncreased levels of residual 
accretion (Fridriksso n et al.l 120101) . Furthermore, long-term 
monitoring of the neutron star transient Cen X-4 in quiescence 
has shown variability in the thermal component that cannot be 
explained by crustal cooling, and may a lso be linked to resid- 
ual accretion during quiescence (ICackett et al.l 12010). The 
monitoring of KS 1731-260 is relatively sparse (only 8 ob- 
servations over 8 years) and thus any increase in flux similar 
to XTE J1701-462 could easily have been missed. However, 
only the first observation of KS 173 1-260 sho ws a significant 
detecti on above 3 keV, and in this observation Wijnan ds et al.1 
(120011) found that the contribution of the power-law was lim- 
ited to about 15% of the flux, though the power-law compo- 
nent was not statistically required to fit the spectrum. In order 
to investigate the potential contribution from power-law emis- 
sion, we created a hardness-intensity diagram using only the 
Chandra count rates in different bands. The data are consis- 
tent with a cooling neutron star atmosphere, though a small 
contribution (a few tens of percent) from a power-law compo- 
nent cannot be excluded. 

From a theoretical perspective, it is puzzling that the de- 
cay is best-fit by a single power law. Numerical simulations 
of crustal cooling expect that the lightcurve should approxi- 
mate a broken power- law which flattens to a constant at late 
times ( IBrown & Cumming 2009). The break in the power- 
law is caused by a change from classical to quantum crys- 
tals, and the timescale on which this occurs is set by the ther- 
mal diffusion time to the depth of this transition. Though in 
KS 1731-260 we only see a single, unb roken power-law de- 
cay in the lightcurve, fitting the model of lBrown & Cum ming 
(2009) allows for quantitative constraints on the thermal con- 
ductivity of the crust. The best-fitting model gives an impurity 
parameter Qi mp = 4.0 and core temperature T c = 5.4 x 10 7 K 
(note that these values are for a fixed mass and radius, and 
any increase in the surface gravity shortens the cooling time 
by decreasing the crust thickness). The impurity parameter is 
marginally higher and the c ore temperature is slightly lower 
than the best-fitting values in lBrown & Cumming (2009f). The 
additional Chandra observation requires a lower core temper- 
ature to fit the lightcurve at late times, where as the impurity 
parameter needs to be increased to maintain the thermal cool- 
ing timescale (a lower core temperature increases the ther- 
mal conductivity and thus increasing the impurity parameter 
compensates). The implied core temperature only requires 
standard models fo r core cooling without enhanced levels 
on ne utrino cooling (Shte rnin et al.l2007tlBrown & Cumming! 
50091) . 

It is interesting to note that a power-law decay with a slope 
of 1/8 (as observed here) is expected for thermal relaxation of 
the core by enhanced neutrino emission (i.e., the rate goes as 
T 6 , see e.g.. iPage et al.ll2006b . Such a scenario requires the 
core to be heated up substantially during an outburst. How- 
ever, in order to do this, or to cool the core over the observed 
timescale, would require a core specific heat several orders of 
magnitude smaller than that provided by degenerate electrons, 
assuming an electron fraction of order 0. 1 . 

EMC thanks the staff at Durham University for kind hos- 
pitality during a visit where much of this work was com- 
pleted. EMC gratefully acknowledges support provided by 
NASA through the Chandra Fellowship Program. 
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